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From Vlasov kinetic equation for collisionless plasmas distribution function in square- 
law approximation on size of electromagnetic field is received. Formulas for calculation 
electric current at any temperature (any degree of degeneration of electronic gas) are 
deduced. The case of small values of the wave numbers is considered. It is shown, that 
the nonlinearity account leads to occurrence the longitudinal electric current directed 
along a wave vector. This longitudinal current orthogonal to known transversal 
classical current, received at the linear analysis. From the kinetic equation with 
Wigner integral for collisionless quantum plasma distribution function is received 
in square-law on vector potential approximation. Formulas for calculation electric 
current at any temperature are deduced. The case of small values of wave number 
is considered. It is shown, that size of a longitudinal current at small values of 
wave number and for classical plasma and for quantum plasma coincide. Graphic 
comparison of dimensionless size of a current quantum and classical plasma is made. 


INTRODUCTION 

Dielectric permeability of quantum plasma was studied by many authors 
m-H- It is one of the major characteristics of plasma also it is applied in 
the diversified questions physicists of plasma 0 - mi 

Let us notice, that in work 111 the formula for calculation of longitudinal 
dielectric permeability into quantum plasma for the first time has been 
deduced. Then the same formula has been deduced and in work [2j. 

In the present work formulas for calculation electric current in classical 
and quantum collisionless plasma at any temperature, at the any degrees of 
degeneration of electronic gas are deduced. 

The approach developed by Klimontovich and Silin [1] is thus generalised. 



At the decision of the kinetic equation we consider as in decomposition of 
distribution function, and in decomposition and sizes of the self-consistent 
electromagnetic held and Wigner integral the sizes proportional to square of 
intensity or potential of an external electromagnetic held. 

Electric current expression consists of two composed. The first composed, 
linear on intensity of an electromagnetic held, is known classical expression of 
an electric current. This electric current is directed along the electromagnetic 
helds. The second composed represents an electric current, which is propor¬ 
tional to an intensity square of electromagnetic held. The second current it is 
perpendicular to the hrst and it is directed along the wave vector. Occurrence 
of the second current comes to light the spent account nonlinear character 
interactions of an electromagnetic held with classical and quantum plasma. 

In works [T5] and [IT] nonlinear effects are studied into plasma. In work 
P] the nonlinear current was used, in particulars, in probability questions 
disintegration processes. We will note, that in work urn it is underlined 
existence nonlinear current along a wave vector (see the formula (2.9) from 

HU)- 


1 CLASSICAL PLASMAS 


Let us show, that in case of the classical plasma described by the Vlasov 
equation, the longitudinal current is generated and we will calculate its 
density. On existence of this current was specihed more half a century ago 
H3- We take the Vlasov equation describing of behaviour of collisionless 
plasmas 

% + bl +e ( E+ c [v ’ H1 )v = 0 ' (L1) 

Electric and magnetic helds are connected with the vector potential by 
equalities 

_ 1 dA iuj . 

E =-— = —A. H = rotA. 

c dt c ’ 

Therefore, 

H = —Ey ■ (0,0,1), [v, H] = —Ey ■ (Vy, ~V X , 0), 

UJ UJ 
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Let us operate with the successive-approximations method, considering as 
small parametre size of intensity of electric held. Let us copy the equation 
(1.1) in the form 

Qf ( k ) Qf( k ) 


dt 


+ V. 


eE v 


uj 


kVy-^- -h (u - kv x ) 

OPx 


X o 

OX 

Qf(k-1) 


dp* 


k = 1,2. 


( 1 . 2 ) 


Here /«» = /„(«) is the absolute Fermi Dirac distribution, 

£ — /m - 1 


/oM = 


1 + exp 


foT J 


mr 


is the electrons energy, y is the chemical potential of electronic 
z 

gas, ks is the Boltzmann constant, T is the plasma temperature. 

Let us consider, that intensity of electric held varies harmoniously: 

E = E 0 e i(kr ' wt) . 

Wave vector we will direct along an axis x: k = k( 1,0,0), and intensity of 

electric held we will direct along an axis y: E = E y { 0,1, 0). 

We notice that _ „ 

<9/o 


v, H 


dp 


= 0, 


because 


<9/o 

dp 


^ V. 


We search for the solution as a hrst approximation in the form 

Z (1) = Zo(P) + Zi. 

where f\ Ey. 

In this approximation the equation (1.2) becomes simpler 


dfi , dh 


dt 


+ u 




= — eE, 


dfo 
' dp y 


(1.3) 
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From (1.3) we receive 


fi = 


2 ieEy P y g(P ) 

Pt u ~ kv T P x ’ 


(1.4) 


where 


n P 2 -a 


g(P) = 


P = P-. 


(l + e p2 -«) 2 ’ p T 

Here vt = a/2 ksT/m is the thermal electrons velocity, px = mvr is the 
thermal electrons momentum, P is the dimensionless momentum. 

In the second approximation for the solution of the equation (1.2) we 
search in the form 

/ (2) = / (1) + h, 


where /\ ~ Ey. 

From the equation (1.2) it is found 


e 2 Ey 

= 2 ~ 

pfpuj 


l r)2 9 ( g(P) 
' VtP v gp x F _ kvi , Pi 


+ 


d(P v g(P )) 

dP„ 


1 


tu — kvxP x 


(1.5) 


Distribution function in the second approximation across the held is con¬ 
structed 

/ = / <2) = / (0) + h + h, (1.6) 

where /i,/ 2 are given by equalities (1.4) and (1.6). 

Let us find electric current density 

2 d 3 p 


j = e vf- 


(1.7) 


(27rh) 3 

From equalities (1.4) - (1.6) it is visible, that the vector of density of a 
current has two nonzero components 

j C jxi jy> 0) • 

Here j y is the density of transversal current, 


Jy 


v vf 


2 d 3 p 

( 27 t ny 


= e 


v vh 


2 d 3 p 
(27rh) 3 


( 1 . 8 ) 
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This current is directed along an electromagnetic field, its density it is 
defined only by the first approximation of function of distribution. The second 
approximation of function of distribution the contribution to current density 
does not bring. 

The density of the transversal current is defined by equality 


ze 


Jy ~ 


2 P 

hjj-i 


Ey(x,t) 


P?g(P)d 3 P 


2iv 3 m y 1 ' J uj — kvrP x 
For density of the longitudinal current according to its definition it is had 




Vxf- 


2 d 3 p 


= e 


= W [ P xh Sp. 


j (27rh) 3 J ‘ (27r/z) 3 (2zr h) 3 

By means of (1.6) from here it is received, that 


9(P y g{P)) , . a / g(^) \ 

dP y T v dP x \ui - kv T P x ' u-korP* 

(1.9) 

In the first integral from (1.9) internal integral on P y it is equal to zero. 
In the second integral from (1.9) internal integral on P x it is calculated in 
parts 


Jx = e 


3 E 2 , 


2 mv^ 


( 27 zhy 


'UJ 


OO 



d / g(P) \ P x dP x 
dP x \uj — kvrP x 2 w — kvrP x 


— OO 



g(P)dP x 
(u - kv T p x y 


Hence, equality (1.9) becomes simpler 

. _ 3 2 2 mvlk f g(P)PyP _ e*E? y q f g(P)PyP 

y (27rh) 3 J (uj — kvrP x ) 3 47 t 3 hmv^ J ( qP x — f2) 3 


Here 


n = 


UJ 


kj^vj 1 


k 

q = ~T 


Equality (1.10) is reduced to one-dimensional integral 


Jx 


OO 


e-'E'iq f ln{l + e a - p *)dP x 


Sir^hmvT 


( q p x - ny 


(i.ii) 


—OO 
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Let us find numerical density (concentration) of particles of the plasma, 
corresponding to Fermi Dirac distribution 


n = / UP) 


oo 

2 dJp 8 ttpt f e a ~ p ~ P 2 dP k^ 




(27rh) 3 (27rh) 3 J 1 + e a p2 2 tt 2 

o 


where Ay- is the thermal wave number, kr = 


ttivt 


OO 

l 0 (a) = J ln(l + e a ~ p )dP. 


In expression before integral from (1.11) we will allocate the plasma (Lang¬ 
muir) frequency 


UJp 


1 4:7re 2 N 


m 


and numerical density (concentration) N 1 and last we will express through 
thermal wave number. We will receive 


J 


■long = E 2 




OO 


q 


ln(l + e a T ~)dr 


y pt 167r/o(a) J (qr — D) 3 


UJp fojJp 

where = -—— =- tt is the dimensionless plasma frequency. 

p k T v T mv% 

This equality we will copy in the form 


jkm S = J classic (0, q)ai t t r kEy 


( 1 . 12 ) 


where aij r is the longitudinal-transversal conductivity, JdasskW, q) is the 
dimensionless part of current, 


T.ir o 


eh ( fojjp \ 2 


-n: 


^classic q) — 


prpKmv^J prkr 

OO „ 

1 f ln(l + e a ~ r )dr 

167r/ 0 (a) J ( qr - D) 3 

— OO 
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The integral from dimensionless part of current is calculated according to 
to known Landau rule 

OO 


ln(l + e a T )dr . 7r 


(qr - ny 


—i 


ln(l + e a_r2 )j " 


— CO 


CO 


+v.p. 


2 q 3 L 

ln(l + e a ~ T2 )dr 
{qr — 12) 3 


+ 


r=Q/q 


—OO 


Symbol V.p. before integral means, that integral it is understood in sense 
of a principal value. 

Let us introduce the transversal electromagnetic field 

„ k(Ek) „ q(Eq) 

&tr — -7"o- ~ & -o-• 

K z q z 

Equality (1.12) can be written down in the invariant form 


j l0ng = Jdassic(^,9)^,trkE^. 


Let us pass to consideration of the case of small values of wave number. 
From expression (1.10) at small values of wave number it is received 


•classic 

Jx 


2e :>, EymVj,k 

(27rh) 3 cj 3 


g(P)Pyd 3 P 


e 3 Eyk%,l 0 (a) 

47T 2 CJ 3 


1 

87T 


e 

mu 



(1.13) 


2 KINETIC EQUATION FOR WIGNER FUCTION 


Shrodinger equation 

ih^- = Hp — H*'p 

dt 

for density matrix p under condition of calibration divA = 0 it will be 
transformed in the kinetic equation [TBJ 




( 2 . 1 ) 
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written down concerning quantum distribution Wigner function 


/(r,P,*)= / P(. r + * pa/fi d 3 a, 


besides 


p( R, R ; , £) = 


1 


/(~T~» P> t)e ip(R_R ' )/fi rf 3 p. 


(2tt^) 3 J j v 2 

Here H is the Hamilton operator, H* is the complex conjugated to H 
operator, H* 1 is the complex conjugated to H operator, acting on the shaded 
spatial variables r'. The scalar potential is considered equal to zero. Integral 
of Wigner is equal (see, example, HU): 

af 


w[f] = 


2 me 


A (r + + A ( r - t) - 2A(r, t) 


dr 


ze 


mch 


a 


a 


A ( r + 2’^) _ A ( r ~ 2 ’^ P ^ + 


ze 


2 mc 2 h - 


A 2 (r + i) - A 2 (r - t)J / } e «(p'-pW»>THT 


(27rh) 3 

Vector potential of an electromagnetic field we take orthogonal to direction 
of the wave vector k (kA = 0) in the form of running harmonious wave 

A(rV) = A 0 e* (kr - Wt) . 

We transform the previous expression of Wigner integral (see, example, 
[II]). We find that 


W[f] = -A(r,t) 


2 me 


_ _ , N _ / hk N N 

V/(r, p - — V) + V/(r, p + — V) - 2V/(r, p, *) 


-A(r,f) 


ze 


mch. 


P 


+A 2 (r,i) ; 


ze 


„/ hk hk ■ 

/(r,P~— ,t)-/(r,P + 

/(r, p - hk, f) - /(r, p + /zk, t) 


2mc 2 h L 

Let us enter local and absolute distributions of Fermi Dirac 


( 2 . 2 ) 


f m = /o(r, C, £) = [! + exp(C 2 


— OL 
























and 


/ (0) = /o(f’) = [l + exp (P 2 ~a)]-'. 


Here 


C = C(r,P,t) = - = P-— A(r,t), 

Vt CPt 


a = 


k B r 


C is the dimensionless electrons velocity, v B = 1/\//^ is the thermal electrons 
velocity, (3 = m/2k B T, P = p/pr is the dimensionless electrons momentum, 
m is the electron mass, k B is the Boltzmann constant, T is the plasma 
temperature, p, is the chemical potential of electronical gas, a is the dimen¬ 
sionless chemical potential. 

Let us show, that the first composed in Wigner integral (2.2) equals 
to zero. We will notice, that according to problem statement gradient of 
quantum distribution function it is proportional to the vector k: V/ ~ k. 
Therefore 


A(r, t) 


/?k Ji k 

V/(r,p - — ,t) + V/(r,p + — ,t) - 2V/(r, p, t) 


~ Ak = 0. 


Thus, Wigner integral is equal 


W[f] = pA 


ze 


ef hk hk . 

/(r, P + —,t) - f(r, p - t) 


-A" 


ze 


2 mc 2 h 


mch . 

2 

/(r, P + hk, t) - /( r, p - /zk, f) 


+ 


Let us return to the kinetic equation (2.1). We will consider convectional 
derivative from this equation 


df 

V dr 



e 


me 


A(r ,t) 


df 


dr 


Thanks to conditions 


— ~ k, Ak = 0 
dr 


we receive that 


v a/ = p_d£ 

dr m dr 




df 
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Let us solve further kinetic Wigner equation for quantum distributions 
function 


df , D d/ , iev T PA 
~di + VTP ~di + ~ch~ 


/( r . p + |a)-/(r,P-|,t) 


-A" 


xe 


2 mc 2 h 


/(r, P + q, t) - /(r, P - ftq, t) 


= 0. 


Here 


hk k , pr 

q = — = 7“, k T = —, 

Pt k‘ r r h 


(2.3) 


kx is the thermal wave number, q is the dimensionless wave number. 


3 SOLUTION OF WIGNER EQUATION 


Let us consider as small parametre size of vector potential of electromag¬ 
netic held A(r ,t). The solution of the equations (2.3) we will be to search 
by the method of consecutive approximations. 

As the first approach for the solution we search in the form, linear on 
vector potential concerning to absolute distribution of Fermi Dirac: 


/ (1) = /o(P) + /i, /i~A(r,t). 


(3.1) 


As the first approximation we take linear concerning to vector potential a 
part of Wigner equation 


<9/ + p df_ + iev T PA 

dt dr mch 


/o(p+f) - /o(p - f: 


= 0. 


Here 


/o(P ± f) = 


q 


m -1 


1 + exp 

Substituting (3.1) in (3.2), we receive the equation 


p ±-l -a 


-i(ui - >’,kPj/, = -^IPA 


ch 


/o(P + f) - /o(P - f) 


From this equation we receive 


, _ m PA / o(p + q/2 ) - /o < p ~ 


ch 


l o — u^kP 


(3.2) 


(3.3) 
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Hence, as a first approximation according to (3.1) and (3.3) the solution 
it is constructed 

/“> = MP) + ^ PA /o(P + q/2)^/o(P^q/2) . 

ch u) — vt kP 

In the second approach we search for the solution in the form 
/ (2) = / (1) + h, /2 ~ A 2 (r, t). 


In the equation (2.3) in the first square bracket function / we will replace 
on /3 1 ', and in the second square bracket function / we will replace on /q, 
i.e. / (2 ^ we search from the equation 


df {2) _^/ (2) iev T PA 

+ n T P—--b 






ch 




q 


-A" 


ze 


2mc 2 h 


/o(P + q)-/o(P-q) 


= 0. 


(3.5) 


In this equation according to (3.4) 


/ (i, (p + fHo(p+f) 

/ (i) ( p -|) =/o ( p -|) 


eww p q\ A /o(P + q) -/o(P) 

ch V 2/ cj — cqT/P + q/2) ’ 

_ q3 A /o(P)-/°( p -q) 

ch V 2/ cj — z>rk(P — q/2) ’ 


We notice that 


Aq = 0, because Aq ~ Ak = 0. 
Therefore the previous two parities become simpler 


/ (1, (P-f)=/o(P-f) 

/ (1) (p + f)=/„(p + f) 


, ^ /QlPl-ZotP-q) 

ch [ uj — nyk(P — q/2) ’ 

, ev T /p A \ /o( p + q ) ~ /o(P) 
ch cj — nrk(P + q/2) 
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The equation (3.5) we will copy in an explicit form 


d/l 1 ) 5/W ievr PA 

+ v T P— -b 


dt dr 

dfo 5/2 

+^ + n T P-^ + 
at or 

/o(P)-/o(P-q) 

cj — nrk(P — q/2) 


ch 


/o(P + f) - /o(P - f ) 


+ 


^4(PA) 2 r/o(P + q) - /o(P) 


c 2 h 2 
ie 2 A 2 


-cj — r>rk(P + q/2) 


/o(P + q) - /o(P - q) 


= 0. 


2mc 2 h 

First three composed in this equation give zero agree to the equation (3.2). 
The rest part of equation leads to equality 


.e 2 4(PA) 2 r/o(P + q)-/o(P) 


\2 r 


— 2i(uj - uykP )_/*2 = — i 
/ 0 (P)-/o(P-q 


cj — nrk(P — q/2). 
From here we find that 


+ 


c 2 h 2 

ze 2 A 2 r 
2 mc 2 h - 


.uj — VT~k(P + q/2) 
/o(P + q)-/o(P-q) 


/2 = 


e 2 n 2 (PA) 2 


7o(P + q)-/o(P) 


2c 2 h 2 (uj — n^kP) Iuj — n T k(P + q/2) 

/o(P) - /o(P - q)l e 2 A 2 / 0 (P + q) - /o(P - q) 


uj — nyk(P — q/2) J 4 mc 2 h uj — n^kP 

So, quantum distribution function of Wigner is constructed and it is 
defined by equalities (3.1), (3.3) and (3.6): / = fo(P) + /1 + / 2 , or, is more 
detailed 

e 2 u 2 (PA) 2 


/ - UP, + + 

ch uj — u^kP 2 c z h z {uj — u^kP) 

;; r /o(P + q)-/o(P) /o(P ~ q) ~ /o(P) l 
Vuj — uyk(P + q/2) n; — uyk(P — q/2). 

e 2 A 2 / 0 (P + q) - / 0 (P - q) 


X 


4 mc 2 h 


uj — u^kP 
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4 ELECTRICAL CURRENT IN QUANTUM PLASMA 


By definition, the electric current density is equal 



(4.1) 


Substituting in equality (4.1) obvious expression for velocity 



and distribution function according to equality / = fo(P) + A + A- 


Leaving linear and square-law expressions concerning vector potential of 
a field, we receive 




(4.2) 


Let us show, that the formula (4.2) for electric current density contains 
two nonzero components: j = (A? AT)- O ne component j y it is linear on 
potential of an electromagnetic field and is directed lengthways field. It is the 
known formula for electric current density, so-called "transversal current". 
The second component j x is quadratic on potential of field also and it is 
directed along the wave vector. It is "longitudinal current". 

The first composed in (4.2) is linear on vector potential expression, and 
second is square-law. We will write out these composed in obvious form 


J 


• linear 



/ S'™ 




(2tt h) 3 



and 


(27rh) 3 
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e 2 ir/P(PA 


2 c 2 h 2 


/q(p + q) — fo{P) 

uj — v T k(P + q/2) 


/o(P)-/ 0 (P-q) l 
uj — , t; T k(P — q/2). 


eVPA 2 

4 mc 2 h 


[/o(P + q) - /o(P - q)] 


PP 

uj — vrkP 


(4.4) 


Expression (4.3) is linear expression of the electric current, found, in 
particular, in our previous work m- This vector expression contains only 
one component, directed along the electromagnetic fields. Really, if a wave 
vector to direct along an axis x i.e. to take k = P(l, 0,0), and potential 
electromagnetic fields to direct along an axis y. i.e. to take 


A(r ,t) = (0, A y (x,t),0), 


from the formula (4.3) we receive 

•linear = 2e 2 p%A y f / fp{P x + q/2) - fp{P x - 
y (27rh) 3 mcq J V P x — Vt/q 


Here 

,/o(Px ± q/2) 


\ _|_ g(A±9/2) 2 +P^ +P^—a 


—Pt + qfo(P))d 3 P. 

(4.5) 


n = 


UJ 

k]'V r j' 


Let us consider expression for an electric current (4.4), proportional to a 
square of potential of an electromagnetic held. Let us notice, that the first 
composed in this expression is equal to zero. Hence, this expression becomes 
simpler 


;quadr 


2 epf, 

(2tt h) 3 


e 2 4P(PA) 2 r / 0 (P + q) -/q(P) 

2c 2 h 2 {uj — n^kP) La; — ryT/P + q/2) 


/o(P-q)-/o(P) l 

uj — nrk(P — q/2). 


e 2 « T PA 2 / 0 (P + q) - /p(P 
4 mc 2 H uj — uykP 



(4.6) 


Let us notice, that vector expression (4.6) contains one nonzero of the 
electric current component, directed along the wave vector 


• quadr _ 
Jx 


e 3 p|A 2 

(2nh) 3 c 2 m 2 VT 


r /o(Px + q) - MP) MP X -q)- /o(P) l 2 
L qP x + q 2 /2-n + qP x -q 2 /2-Q J y + 
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(4.7) 


[MP* + q)~ ,/i'i( P, - (/)] 

Let us lead to a kind convenient for calculations, the formula (4.7) for 
density of a longitudinal current. 

Let’s consider the first integral from (4.7). We will calculate the internal 
integrals in plane ( P y , P z ), passing to polar coordinates 

OO 

I f 0 (P x ± q)PydP y dP s = wf pln(l + e^^-^dp, 

0 


P x d 3 P 
qP x ~ A 


J fa{P x ± q)dP y dP = = n ln(l + e^^). 

Thus, size of density of generated longitudinal current into quantum plasma 
it is equal 


•quant _ 
Jx 


Tie 3 p 3 T A 2 
(2n h) 3 m 2 c 2 VT 


OO 



—OO 


L{P X + q, P x ) L(P X — q, P x ) 
qP x + q 2 /2 - n qP x - q 2 /2 - U 


1 \ _|_ e a-(P x +q) 2 

H-In - — -Tn 

2 1 + e a-{P x -q) 2 


P x dP x 


Here 


OO 


L(P X db q, P x ) = J pin 


qP x ~ A 

l _)_ g a-(P x ±q) 2 -p 2 
l _|_ e a-P2-p 2 


(4.8) 


dp. 


Let us transform the expression standing in integral (4.8). At first let us 
pass from potential to intensity of a held A y = —( ic/uj)E y . We will receive 


7re 


3 p 3 T A 2 


7re 3 k^E 2 


{2tt h) 3 m 2 c 2 VT 87 PnPvru 2 


Let us transform this expression by means of expression for the thermal 
wave numbers W- We receive, that 

7 ze 3 k^E 2 e 3 NE 2 1 e /cj p \ 2 ^ 2 

8jr 3 m 2 VTU 2 4:lo(a)m 2 VTU 2 167T itivt V c u J y 
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1 


1 e /Q p \ 2 2 _ 
167TPT \ Q j V 


1 


V -El = 


167rlo(a)qQ, 2 l ' tr v ‘ 


167 r/o(ct)f 2 2 pt y 

where the quantity of longitudinal-transversal coductivity oy ir was introduced 
in part 1: cr Z)tr = eV? p /{p T k T ). 

Now equality (4.8) we will present in the form 

Jx Uailt = ^quant(^, q)&l,trkEy(x, t). (4.9) 

In (4.9) Jq U ant(^ 5 g) is the density of dimensionless longitudinal current, 

OC 

1 f 


'/quant(^4? 0 ) — 


167r/o(ct)qkL 


1 1 + e a -( r+ «) 

H— In- 7 -r 

2 1 + e a -( T "9) 


L(r + q : r) L(t - q,r) 
qr + q 2 /2 — Q qr — q 2 /2 — 

rdr 


qr — Vt 

Let us transform integral from (4.10) from the first composed 


(4.10) 


oo r 


Ji = 


L{r + q,r) L(r — q,r) 


qr J rq 2 /2 — Q qr — q 2 /2 — Vt 


— OO 


rdr 
qr — £1 


= J 2 + J‘i- 


Here 


OO 


J2 = 


L(r + q, r)rdr 


—OO 
OO 


Js = 


(qr + q 2 /2 - Q) (qr - Q) ’ 

L(r — q, r)rdr 
(qr — q 2 /2 — Q)(qr — Q) 


—OO 


In integral J 2 we will make variable replacement r —>■ r — qf 2, and in 
integral J 3 we will make replacement r —> r -\- q/2. As result it is received 


J 2 = 


^3 = 


L(r + q/2,r — q/2)(r - q/2)dr 
(qr — 0,)(qr — q 2 /2 — Q) 

O 

f L{r - q/2, r + g/2)(r + q/2) 
(qr - f 2 )(qr + q 2 /2- Q) 


—OO 
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We notice that 


7 1 -)_ e a-(r-g/2) 2 -p 2 

L{r-q/2,T + q/2) = / pin „ , n _ UMa/0 w_^ dp = 


OO 

= -J P In 
0 


X _)_ e a-(r+g/2) 2 -p 2 

0 

X -1- e a-(r+ 72 ) 2 -p 2 

—jhp = —L(t + g/2, r - q/2). 


X -|- e a-(T-q/2) 2 -p 

Hence, integral J\ equals 

L(r + q/2,T - q/2) r r - q/2 


OO 


7i = 


—OO 


gr — ft 


= (/ft 


r+ g/2 


qr — q 2 / 2 —ft qr + q 2 / 2 —ft 


L(t + q/ 2, r — q/2)dr 


dr 


(qr - n)[(qr - ft) 2 - q A / 4]' 
Let us calculate the second integral from (4.10) 


— OO 


OO 


1 f 1 + e“-( r+ 9) rdr 
4 ~2J ln 1 + e «-(r-W qT -n 

—OO 

= 1 _JL_ ldT = 

2q J 1 + e a- ( r- 7 V qr — Q. 

— OO 

00 / x9 

ft r i + e a -( T+,? ) hr 

~ 2g J n 1 + e a_ ( r_ 7 2 gr - ft 

—OO 


ft 


OO 


= 2g J ^ ) 


1 


1 


.qx — q 2 — ft qx + (/ 2 — ft. 


dx = 


— OO 


_ f ln(l + e a x2 )hr , 

= ft(/ / t-—-- —dx. 


{qx — ft) 2 — g 4 /4 


Finally, integral (4.10) equals 

J q ua„.(a?) = - 16 ^ o(a)f}! 


7(r + g/2, r - q/2)dr 
(gr-ft)[(gr-ft) 2 -g 4 /4] 


— OO 
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+ 


ln(l + e a x2 )dx 
(qx — 12) 2 — q 4 /4 


dx 


—OO 


(4.11) 


At calculation singular integral from (4.10), which not writing out let us 
designate through 7(12, q), it is necessary to take advantage known Landau 
rule. Then 

7(12, q) = Re 7(£1, q) + i Im 7(12, q). 

Here 

OO 

L(t + q/2, t - q/2)dr 


Re 7(12, q) = V.p. J 


—OO 
OO 


+v.p. 


(gr - 12) [(gr - 12) 2 - g 4 /4] 
ln(l + e a ~ x2 )dx 


dx , 


— OO 


(qx — 12) 2 — g 4 /4 
Symbol V.p. means principal value of integral also, 

Im 7(12, g) = (q 2 - 212)l (^ + ^) + 


<r 


q 


+(g 2 + 212)L(-| + ^,| + ^) + 


+212 




r . 12 12\ /12 12 + 12</ 2 , l + e «-(V9+?) 2 

L I ~q + iJ " L W ” '^ J + ~T '1 + e-CV,-,)’ 


Considering antisymmetry of function L(X, Y ) on the variables (L(X, Y ) = 
— L(Y : V)), we will simplify expression for an imaginary part of density of 
the dimensionless current. It is as a result received, that 


Im 7(12, q) = 


7t12 

q 4 


ATI 'q V q 12 \ 

-4L(T + -,-^ + - + 
-2 q 2 qJ 


+212 


+ /12 12 \ /12 12 \ 


12o 2 1 + e U-(n/q+qY 

H-—In 


1 _p e a-{Q,/q-qf 


q qJ. 

Equality (4.10) for density of a longitudinal current we will present into 
the vector form 

j quant = J quanl (n j 9 ) CTiifr feE 2 .. 


18 

















Let us show, that at small values of wave number density longitudinal 
current both in quantum and in classical plasma coincide. 

According to (4.7) at small q after linearization fo(P x ± q) it is received 

, 2 e 3 p?rAlq f 0 „ 27re 3 jAA 2 /o(ad 

■quant _ _ £_£ U 1 / p\p2rf3 p _ A 

Jx (2tt h) 3 c 2 m 2 u T 0 J J x {2irhfc 2 m 2 u 


Let us take advantage of relation of thermal wave number with the nume¬ 
rical density, and also relation of potential and intensity of electromagnetic 
held. We receive, that 


•quant 

Jx 


1 

87T 


e 


muj 



(4.11) 


Expression (4.11) in accuracy coincides with (1.12). These expressions let 
us copy in the vector kind 



87t muj \uj ) 


5 CONCLUTIONS 

Graphic investigations of sizes of density of longitudinal current we will 
spend for case: chemical potential equals to zero (a = 0). Curves 1 answer 
to classical plasma, curves 2 to quantum. 

On figs. 1 and 2 we will present behaviour real (fig. 1) and imaginary 
(fig. 2) parts of density of dimensionless longitudinal current at Q = 0.5 
Depending on dimensionless wave number q. At decrease parametre Q curves 
1 and 2 approach and become indiscernible. 

On fig. 3 and 4 we will represent behaviour real (fig. 3) and imaginary (fig. 
4) density parts of longitudinal current depending on dimensionless wave 
number q in the case 0 = 1. From these drawings it is visible, that with 
growth dimensionless frequency of oscillations of an electromagnetic held of 
ordinate graphics quickly decrease. 

On hg. 5 and 6 we will represent behaviour real (hg. 5) and imaginary 
(hg. 6) parts of longitudinal current depending on dimensionless frequency 
of oscillations of electromagnetic held 0 in the case q = 0.3. At reduction 
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dimensionless wave number q curves 1 and 2 approximate and at small q 
practically coincide. 

On figs. 7-10 we sonsider behaviour of classical plasmas. 

On fig. 7 and 8 the behaviour real (fig. 7) and imaginary (fig. 8) is 
represented parts of longitudinal current depending on dimensionless wave 
number, 0. Curves 1,2,3 accordingly answer values dimensionless chemical 
potential a = —5, 0, 3. 

On fig. 9 and 10 we will represent behaviour real (fig. 9) and imaginary (fig. 
10) parts of longitudinal current depending on dimensionless wave number q 
in a case 0 = 1. Curves 1,2,3 answer according to values of the dimensionless 
chemical potential a = —5, 0, 3. 

In the present work the account of nonlinear character of interaction 
electromagnetic field with classical and quantum plasma is considered. It has 
appeared, that the account of nonlinearity of an electromagnetic field finds 
out generating of an electric current, orthogonal to a direction of electromag¬ 
netic field. 

Further authors purpose to consider new problems about oscillations plas¬ 
mas and about skin-effect with use square-law onto potential of expansion of 
distribution function. 
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Fig. 1. Real part of density of dimensionless current, Id = 0.5. 



Fig. 2. Imaginary part of density of dimensionless current, fl = 0.5. 
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Fig. 3. Real part of density of dimensionless current, 0 = 1. 



Fig. 4. Imaginary part of density of dimensionless current, 0 = 1. 
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Fig. 5. Real part of density of dimensionless current, q = 0.3. 



Fig. 6. Imaginary part of density of dimensionless current, q = 0.3. 
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Fig. 7. Real part of density of dimensionless current in classical plasma, q — 1. Curves 
1,2,3 correspond to values of dimensionless chemical potential a = —5, 0, 3. 



Fig. 8. Imaginary part of density of dimensionless current in classical plasma, q — 1. 
Curves 1,2,3 correspond to values of dimensionless chemical potential a = —5,0,3. 
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Fig. 9. Real part of density of dimensionless current in classical plasma, 0 = 1. Curves 
1,2,3 correspond to values of dimensionless chemical potential a = —5, 0, 3. 



Fig. 10. Imaginary part of density of dimensionless current in classical plasma, 0 = 1. 
Curves 1,2,3 correspond to values of dimensionless chemical potential a = —5, 0, 3. 
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